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r ^ ' Abstract. We consider four-dimensional Riemannian manifolds with commuting higher 

P^H , order Jacobi operators defined on two-dimensional orthogonal subspaces (polygons) and 

• on their orthogonal subspaces. 

^v^j ' More precisely, we discuss higher order Jacobi operator J{X) and its commuting 

associated operator J{X • ) induced by the orthogonal complement X of the vector X, 
i. e. J(X) o J{X X ) = J{X X ) o J(X). 
■ At the end some new central theorems have been cited. The latter are due to P. 

' Gilkey, E. Puffini and V. Videv, and have been recently obtained. 

•5 

1. Preliminaries 

Let (M, g) be a n-dimensional Riemannian manifold with a metric tensor g. Tangent 
£f) ■ space at a point p £ M we denote by M p , and let S p M be the set of unit vectors in M p , 

£> | i. e. S P (M) := {z £ M p \ \\g(z,z)\\ = 1}. Let T(M) be the algebra of all smooth functions 

on M and X{M) be the JF(M)-module of all smooth vector fields over M . Let also 

R(X, Y)Z = VxVyZ - VyVxZ - V [x ,y]Z, X,Y,Z£ X(M) 

' be the (1, 3) curvature tensor of the Levi-Civita connection V. We define 



o 



R(X, Y, Z, U) := g(R(X, Y, Z), U) 
to be the associated (0, 4)-curvature tensor which satisfied the following algebraic properties: 

i) R(X,Y,Z,U) = -R{Y,X,Z,U), 

ii) R(X, Y, Z, U) = -R(X, Y, U, Z), 

iii) R(X, Y, Z, U) + R(Y, Z, X,U)+R(Z, X, Y,U)=0 (first Bianchi identity) , 

iv) R(X,Y,Z,U) = R(Z,U,X,Y). 

In the Riemannian geometry the following differential equality is also true: 
v) (V X R){Y,Z,W) + (V Y R)(Z,X,W) + {V Z R)(X,Y, W) = (second Bianchi identity), 
where 

{VxR)(Y, Z, W) := V X (R(Y, Z)W) - R{V X Y, Z)W - R(Y, V X Z)W - R(Y, Z)V X W, 

and V X R is the covariant derivative of the (1, 3)-curvature tensor R with respect to X, 
X,Y,Z £ X{M). 
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Let J : M p — > M p be the Jacobi operator defined by: 

(1.1) J(X)U = R(U,X,X). 

One can easily see that J(X)X = and g(J(X)Y, Z) = g(Y,J(X)Z) which means that 
Jacobi operator is a symmetric linear operator. 

Jacobi operator can be diagonalized in the Riemannian geometry. In this case we say that 
g is Osserman metric if the eigenvalues of the Jacobi operator are constant over the tangent 
bundle S(M) := {J peM S p M. If (M, g) is a rank one locally symmetric space, i. e. VR = 0, 
where V is the connection with all positive or all negative sectional curvatures [2] or (M, g) 
is flat, i. e. R = 0, the group of local isometries acts transitively on S(M) and each Jacobi 
operator has constant eigenvalues. Osserman [5] conjectured that the opposite is also true 
and this was confirmed by Chi when n = 4 and n = (2 mod 4) [2] and by Nikolaevsky [7] 
when n ^ 16. 

Gilkey, Stanilov and Videv [5] introduced a new operator which they called general Jacobi 
operator of order k or k-order Jacobi operator. More precisely, if {1^}*^ is any orthonormal 
basis for an arbitrary fc-plane ir € M p , the higher /c-order Jacobi operator is defined by: 

(1.2) J(tt):Y^ ]T R(Y,Yi)Y i= £ J(Y t )Y. 

l<i<k l<i<k 

It can be easily verified that this operator does not depend on the basis of 7T. 

2. Some commutativity conditions 

Another variety of problems, connected to the higher order Jacobi operator, emerged 
thanks to Stanilov and Videv [11] . They are connected with some commutativity conditions 
forced on (1.2). Recently Brozos- Vazquez and Gilkey [1] were able to prove the following 

Theorem 2.1. Let (M,g) be a Riemannian manifold, dimM > 3. Then 

(A) (M,g) is flat iff J{X)J(Y) = J{Y)J{X) for arbitrary vectors X, Y <E M p ; 

(B) (M, g) is a manifold with a constant sectional curvature iff J(X)J{Y) = J(Y)J(X) 
for arbitrary vectors X, Y £ M p such that X JlY . 

In this paper authors will characterize indecomposible four-dimensional Riemannian man- 
ifolds that satisfy the following two conditions: 

For arbitrary unit vector X £ M p , p € M , we have: 
(CI) J{X) o J{X ) = J{X i -) o J(X), where X 1 - is the orthogonal complement of X in 
M p . 

For arbitrary 2-plane a C M p , p 6 M , we have 
(C2) J{a) o Jfa 1 ) — J(a^)oJ(a), where is the orthogonal complement of a in M p . 
Our main goal is to prove the following 

Theorem 2.2. Let (M,g) be a four- dimensional indecomposible Riemannian manifold. 
Then the following are equivalent: 

(a) Equality (CI) holds for arbitrary unit vector X € M p , p € M ; 

(b) Equality (C2) holds for arbitrary 2-plane a C M p , p € M p ; 

(c) (M, g) is Einstein. 
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Proof, (a) => (c) Let {ei, ei, e3, e4} be an ortho-normal basis for M p , p e M. Curvature 
operator matrices i7{ ei , e2 , e3 } and J7{ e4 } then have the form: 



(2.1) 



and 



(2.2) 



where p l j = p(e,, ej 



/K 12 + K 13 




#1332 


#1223 




P14 


\ 


#1332 


if 12 + ^ 13 


#2113 




P2A 




#1223 




#2113 


#13 + if 23 




P34 




V P14 




P24 


P34 


if 14 + if 24 


+ if 34 / 






/ #14 


#1442 #1443 












#1442 


if 24 #2443 













#1443 


#2443 if 34 
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y~^ff(i?(efc, ej)ej,ek) are the components of the Ricci (0,2)-tensor 
1,...,4. 



fc=i 



p and ify := g(R(e l , ej)ej, e 4 ), i, j 
We have the matrix equality 

( 2 - 3 ) ^{ei,e 2 ,e 3 } ^{e 4 } = ^{e 4 } ^{ei,e 2 ,e 3 }> 

which leads us to the equations: 

(ei) #14(#1224 + #1334) + #1442(#2114 + #2334) + #1443 (#3114 

(e 2 ) #1442 (#1224 + #1334) + #24(#2114 + #2334) + #2443 (#3114 
(63) #1443 (#1224 + i?1334) + #2443 (#2114 + #2334) + #34(#3114 



#3224) 
#3224) 
#3224) 



We do a cyclic change of indecies 1 
the equations 



1 in (ei), (e 2 ), (e 3 ), and get 



(ei) ifl2(i?1332 + #1442) + #2113(#1223 + #1443) + #2114 (#1224 + #1334) = 0, 
(el) #2113(#1332 + #1442) + ifl3(#1223 + #1443) + #3114(#1224 + #1334) = 0, 
(e\) #2114(#1332 + #1442) + #3114(#1223 + #1443) + #u(#1224 + #1334) = 0. 



Doing a cyclic change of indecies 1 



1 in (ej), (e?) and (el), we get 



(e\) if23(#2113 + #2443) + #3224(#2114 + #2334) + #1223(#1332 + #1442) = 0, 
(e 2 ) #3224(#2113 + #2443) + if24(#2114 + #2334) + #1224(#1332 + #1442) = 0, 
(e|) #1223(#2113 + #2443) + #1224(#2114 + #2334) + ifl2(#1332 + #1442) = 0. 



Another cycling change 1 



1 in (el), (e|) and (e|) will give us 



( e l) ^34 (#3114 + #3224) + #1334(#1223 + #1443) + #2334(#2113 + #2443) = 0, 
( e 3) #1334(#3114 + #3224) + ifl3(#1223 + #1443) + #1332(#2113 + #2443) = 0, 
(e 3 ) #2334 (#3114 + #3224) + #1332 (#1223 + #1443) + if 23 (#2113 + #2443) = 0. 



Solving (ei), (e 2 ), (e 3 ) with respect to #1224 + #1334, #2114 
Maple, we get: 



12113 T -T12443J 
#2334, #3114 " 



#3224, using 



#1224 + #1334 — #2114 + #2334 — #3114 + #3224 — 0, 
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since the above is in fact the trivial solution to the system of equations (ei), (e2), (e3) which 
is homogeneous: 

#12 R\ii2 R\AA3\ I -Rl224 + i?1334\ / 0\ 
i?1442 # 24 i?2443 I I ^2114 + i?2334 I = I I . 
i?1443 -^2443 #34 / \-R3ii4 + #3224/ \0/ 

Analogously, solving the other two homogeneous systems (e\), (el), (ef) and (eg), (e§), 
(ef) we get that 

i?1332 + i?1442 = i?1223 + -Rl443 = R\22i + i?1334 = 0, 

and 

-Rl223 + i?1443 = i?2113 + i?2443 = #3114 + i?3224 = 0. 

From (2.3) we also get 

(ei) (if 14 — #24j#1332 + -^2113-^1443 — i?2443-Rl223 = 0, 
( e s) (#14 — #34)-Rl223 + i?2113-Rl442 — i?2443#332 = 0, 

( e 6) (#24 — #34)i?2113 + (#13 — if I2)i?2443 + -^1223-^1442 — i?1332-Rl443 = 0. 

By a cycling change of indecies 1 — > 2 — > 3 — ► 4 — > 1 in (ei), (es) and (e^), we get 

(e\) (if 12 — ifl3)i?2443 + -^2114-^3224 — #3114^2334 = 0, 
(el) (if 12 — ifl4)i?2334 + i?3224i?2113 — -^3114-^2443 = 0, 

(e 4 ) (if 13 — if I4)i?3224 + (if 24 — #23)#3114 + i?2334#2113 — i?2443-R2114 = 0. 



Repeating the same procedure two more times, we get 



(el) 


(#23 


— #24)#3114 " 


- -Rl334-Rl223 


— #224#443 = 0, 






(el) 


(#23 


— #12)-Rl443 " 


~ i?3224i?1334 


— i?1224#3114 = 0, 






(el) 


(if 24 


— #12)i?1334 " 


1- (#13 - #34 


)i?1224 + i?1443-R3224 


— i?3114i?1223 


= 


and 














(el) 


(#34 


— #13)^1224 - 


" i?2334-Rl442 


— #332-^2114 = 0, 






(el) 


(#34 


— if 23)^2114 " 


" -Rl334-Rl442 


— #332-#224 = 0, 






(eg) 


(#13 


— if 23j#1442 - 


1- (#24 - #14 


)i?1332 + i?1334i?2114 


— i?2334i?1224 


= 



Further, from (e 4 ), (e 5 ), (e 6 ); (ei), (el), (el); (el), (ef), (ef); (ef), (ef), (ef) and using 
that 



#2113 - 


- i?2443 


= 0, 


#332 " 


f #442 


= 


-R1223 - 


" #443 


= 0, 


-#224 " 


" #334 


= 


-R2114 - 


" i?2334 


-0, 


#3114 " 


" #3224 


= 0, 



we get the system of equations 



#12 


-#34 


#13 


= #24 


#14 


= #23 



with respect to the basis {ei, e2, e3, 64}. The latter is equivalent to (M, g) being an Einstein|10j.D 
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(c) =>• (a) Suppose (M, g) is a four-dimensional Einstein manifold and let X 6 M p , 
p G M and X 1 - is the orthogonal complement of X. Then p = Aid, A = const., and hence 

J(X) o J(X^) - J(X^) o J(X) = J{X) o J(X X ) + J(X) o J(X) - J{X) o J{X) - J(X^) o 
J{X) o [J(X X ) + J(X)] - [J(X) + J(X X )] oJ(X) = J(X) op-po J(X) = 

A( l 7(X)oId-IdoJ(X)) = 0. 



Analogously, one can prove, using [9], the following 

Corollary 2.1. Lei {M,g) be a three-dimensional Riemannian manifold. Then the next two 
conditions are equivalent: 

(i) J(X) o J{X^) = J{X X ) o J(X) for arbitrary X, X 1 - G M p , p G M. 

(ii) (M,g) has a constant sectional curvature k such that R(X,Y, Z) = n{g(Y 1 Z)X — 
g(X,Z)Y), X,Y,Z eM p . 

(b) (c) Let {ei, e 2 , e^, e A } be an orthonormal basis for M p , p G M. Then the curvature 
operator matrices J^{ ei .e 2 } an d ^{e 3 ,e 4 } have the form: 



(2.4) 



<J{ei,e 2 } ~ 



fK 12 + K 13 


-R1223 






-R1223 


if 12 


-R2113 


-R2113 


if 13 -1- if 23 


-R2114 


P34 



i?1224 
i?2114 
P34 



and 



(2.5) 



(K 13 

P12 
-R1443 

\ i?1334 



if 14 



P12 


i?1443 


if 23 + if 24 


i?2443 


i?2443 


if 34 


i?2334 






R 



1334 
i?23 


#34 J 



Let also 



/if 



12 



(A) 



+ ifl 3 





i?1223 





if 12 


i?2113 


1223 


i?2113 


if 13 + if 23 


1224 


i?2114 


P34 



i?1224 
i?2114 
P34 

K 1A + K 2A J 



\ fK 13 + K u 
P12 

i?1443 
i?1334 



V 



P12 

if 23 + if 24 
i?2443 
i?2334 



i?1443 i?1334\ 

i?2443 i?2334 
if 3 4 

if 34 / 
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Simple computations for the matrix (A) give us: 



On = Ki 2 (Ki 3 + if 14) + i?i223-Rl443 + -^1334^1224, 
«12 = if 12^12 + #1223^2443 + ^2334^1224, 
dl3 = ifl2-Rl443 + K34R1223, 
<X\4 = if 12^1334 + ^34-^1224, 
«21 = ifl2Pl2 + ^1443-^2113 + ^2114-^1334, 
a 2 2 = if 12 (if 23 + if 24) + ^2443-^2113 + -^2114^2334, 
023 = if 12-^2443 + if 34-^21 14, 
a 2 4 = if 12^2334 + -^34-^2114, 
(2.6) a 3 i = (if 13 + K 13 )R 12 23 + (if 13 + if23)-Rl443 + P12-R2113 + /?34-Rl443, 

a 32 = (if 23 + if24)i?2113 + (if 13 + if 23)^2443 + Pl2-Rl223 + ^34-^2334, 
a 33 = (if 13 + if 23)^34 + -R1223-R1443 + -^2443-^2113, 



a 34 = P34if34 + i?1223-Rl334 + -^2334-^2113, 

a 4 i = (if 13 + ifi4)i?1224 + (if 14 + if24)i?1334 + £12-^2114 + P34-TC 
a 42 = (if 23 + if24)i?2114 + (i?14 + if24)i?2334 + Pl2-Rl224 + P34-R: 
a 43 = P34K34 + i?1443-Rl224 + i?2114i?2334, 
a 44 = (if 14 + if 2 4)if34 + i?1334-Rl224 + i?2114i?2334 



-2114 + P34-Rl443, 
•2443) 



On the other hand let 



(B) 



/if 13 + if 14 
Pl2 
i?1443 
\ i?1334 



Pl2 
if 23 + if 24 
i?2443 
i?2334 



i?1443 
i?2443 
if 34 





^?1334\ 
i?2334 


if 3 4 J 



( if 12 


i?1223 
\i?1224 





if 12 

i?2113 
i?2114 



i?1223 
i?2113 
if 13 + if 23 
P34 



i?1224 \ 
i?2114 
P34 



For the matrix (B) elements we derive: 
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(2.7) 



bll = Ki 2 (Ki 3 + K14) + -Ri223-Rl443 + ^1334-^1224, 
bl2 = K12P12 + i?1443-R2113 + -^2114-^1334, 

^13 = (-^13 + -^14)^1223 + (-^13 + -^23)^1443 + P12-R2113 + P34-Rl334, 
&14 = (-^13 + ^14)-Rl224 + {K14 + i^24)-Rl334 + ^12-^2114 + P34-Rl443, 
^21 = K12P12 + i?1223-R2443 + ^2334^1224, 

622 = ^12(^23 + -^24) + ^2443-^2113 + ^2114-^2334, 

623 = (-^23 + K2A)R2Wi + {K\3 + -^23)^2443 + /0l2-Rl223 + ^34-^2334, 

624 = (-^23 + -f^24)-R2114 + (-^14 + -f^24)-R2334 + Pl2-Rl224 + A34-R2443, 
^31 = -^12-^1443 + K 34 Ri223, 

632 = -^12-^2443 + ^34-^2113, 

633 = (-^13 + K23)K 3 4 + -Ri223-Rl443 + ^2443-^2113, 
^34 = P34-K34 + -Rl443-Rl224 + ^2114-^2443, 

641 = -^12-^1334 + -?^34-Rl224, 

642 = -^12-^2334 + -^34-^2114, 

^43 = P34-^34 + -^1223-^1334 + ^2334^2113, 

644 = (K14 + K 2 i)K si + -R1334-R1224 + ^2114^2334- 

From [C2] we have a\2 = &12 and according to (2.6) and (2.7), we get: 

(ej) i?1223-R2443 + ^2334-^1224 _ ^1443^2113 — -^2114-^1334 = 0, 
(es) -Rl223-Rl334 + ^2334-^2113 — ^1443^1224 — -^2114^2443 = 0. 



By a cycling change of indecies 1 



1 in (e?) and (es) , we get 



{ e \) -^2334-^3114 + ^1443-^1332 — ^2114^3224 — -^1223^1442 — 0, 

(eg) -R2334-Rl442 + -^1443-^3224 — -R2114-Rl332 — -Rl223^3114 = 0. 

We solve, using Maple, (e?), (e 8 ), (e\) and (e|) together and arrive at the homogeneous 
system 



(2.8) 



/ -R2443 ^1224 —^2113 

-Rl334 -^2113 —R\22i 

-i?1442 Rzili -Rl332 

\ — -^3114 -Rl442 -^3224 



--Rl334 
i?2443 
^3224 



\ 


^1223^ 




( \ 




-^2334 









-Rl443 







I 


\-R2114/ 




w 



Solving (2.8) with respect to -R1223, -R2334, -R1443 and -R2114, we get at least the trivial 
solution: 



(2.9) -R2114 — -R1223 — -R1443 — -R2334 — 0. 

We also have J{ ei ,e 3 } = ^{e 2 ,e 4 } and J{ ei ,e 4 } = ^{e 2 ,e 3 }> and using (2.9), it follows that 
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(2.10) 



#2113(#1332 + #1442) + #1334(#3114 + -R3224) — 
#1334 (#12 - if 14 - if 24) + i?i 2 2 4 (#34 - if 13 - #14) = 
#2443 (#12 — if 13 — #23) + #2113 (#34 — # 23 — if 24) = 
#1224(#1332 + #1442) + #2443 (#3114 + #3224) = 

-^3224(^2113 + #2443) + #1442 (#1224 + #1334) = 

#144 2 (#23 ~ 12 - if 13) + #133 2 (#14 ~ if 34 - if i3 ) = 
#1332(#2113 + #2443) + #3114(#1224 + #1334) = 

We solve (2.10) with respect to the tensor R components #1332, #1442, #3114, #3224, 
#1224, #1334, #2113, #2114 using, for example Maple, and as a result we get 

i?1224( — if 34 + if 13 + if 14) = #1334(#12 — if 14 — if 24) 
#2113 (-#34 + if 23 + if 24) = #2443 (#12 - if 13 ~ if 23) 

and 



(2.11) 



(2.14) 



(2.12) #1332 = #3224 = #1442 = #3114 — 0. 

Further, by changing the orthonormal basis {ei 1 e 2 ,e 3 ,e^} with the orthonormal basis 

ei-e 2 ei+e 2 e 3 - e 4 e 3 + e 4 1 . , , , . . 
= — , = — , = — , = — > and using 2.11) and (2.12) we receive a new system 

y/2 y/2 V2 y/2 J 

of equations with respect to the tensor curvature components which is equivalent to (2.10). 

From there we can conclude that 

(2.13) i?i334 = #2443 = 0. 

From (2.9), (2.11) and (2.12) it follows that all of the components Rijjk are equal to zero 
for all i, j, k =1,2, 3, 4. 

If we reformulate the second and third equation in (2.10) by changing the basis as shown 
above and transforming the curvature components using (2.9), (2.11) and (2.12), we get 

(if 13 - if 23 + if 14 - #24)(#1432 + #1342 + ifl2) = 
(if 13 + if 23 - if 14 - if 24) (#1432 + #1342 + if 34 ) = ' 

By a cycling change of indecies 1 — > 2 — > 3 — > 4 — > 1 above and some extra tedious 
computations, we get the system: 

(if 13 + if 14 -if 23 -if 24) (if 12 
(if 12 + if 14 - if 23 - if 34) (if 13 
(if 12 + if 13 -if 24 -if 34) (if 14 
(if 12 + if 24 - if 13 - if 34) (if 23 
(if 12+ if 23 -if 14 -if 34) (if 24 
(if 13 + if 23 "if 14 -if 24) (if 34 

Solving (2.14) and (2.15) with respect to the sectional curvature components if 12, if 13 
and if 14, it follows that 

(2.16) #u = #23, if 13 = if 24, if 12 = if 34, 

and since the basis {ei, e 2 , e 3 , e^} has been arbitrary chosen in M p , it follows that (M, g) is 
Einstein [10]. □ 



(2.15) 



#1342 — #1423) = 
#1234 — #1423) = 
#1234 — #1342) = 
#1234 — #1342) = 
#1234 — #1423) = 
#1342 — #1432) = 
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(b) => (a) If (M,g) is an Einstein manifold then p — Aid, A = const., and if a is a 
2-plane in M p , p £ M, it follows that 

J(a) o J(a x ) - J(a^) o 7(a) = 7(a) o J{a^) + 7(a) o J(a) - 7(a) o 7(a) - J{a^) o J{a) = 
7(a) o [J{a x ) + J(a)] - [j(a) + 7(a x )] oJ(a) = 7(a) op-po J(a) = 

" v ' " v ' 

p p 

X (J (a) old -Id 0.7(a)) =0. 

That completes the proof. □ 

3. New approaches and results 

Recently Gilkey, Puffini and Videv [3] were able to generalize the results above. They 
define 971 := (V, (•, •), A) to be a 0-model if (•, •) is a non-degenerate inner product of signature 
(p, q) on a finite dimensional vector space V of dimension m = p + q and if A £ ® A V* is an 
algebraic curvature tensor. 

Let Gr r s (y, (•, •)) be the Grassmannian of all non-degenerate linear subspaces of V which 
have signature (r, s); the pair (r, s) is said to be admissible if and only if Gr TtS (V, (■, •)) is non- 
empty and does not consist of a single point or, equivalently, if the inequalities < r < p, 
< s < q, and l<r + s<m — 1 are satisfied. Let [A, B] :— AB — BA denote the 
commutator of two linear maps. Then they establish the following result: 

Theorem 3.1. Let 971 = (V, (•, ■), A) be a 0-model. The following assertions are equivalent; 
if any is satisfied, then we shall say that 9JI is a Puffini- Videv 0-model. 

(1) There exists (ro,so) admissible so that 
JWoJt^-Jfi^oJW for all tt £ Gr ro , So (V, (;■)). 

(2) 7(7r) o 7(tt ± ) = ,7(71) 7(7r) for every non- degenerate subspace tt. 

(3) [7(tt),p] = for every non-degenerate subspace tt. 

We say that 9Jt = (V, (■,■), A) is decomposible if there exists a non-trivial orthogonal 
decomposition V = V\ ® V2 which decomposes A — A\ ® A2 ; in this setting, we shall write 
9JI = SDti 97l 2 where 971^ := (Vi, (•, -)|yj , Ai). One says that OJl is indecomposible if Wl is not 
decomposible. 

By Theorem 3.1, any Einstein 0-model is Puffini- Vide v. More generally, the direct sum of 
Einstein Puffini- Videv models is again Puffini- Videv; the converse holds in the Riemannian 
setting: 

Theorem 3.2. Let 371 = (V, (•, ■), A) be a Riemannian 0-model. Then Wl is Puffini-Videv 
if and only if 971 = 97ti • • • © 97tfc where the 971^ are Einstein. 

In the pseudo- Riemannian setting, a somewhat weaker result can be established. One says 
that a 0-model is pseudo-Einstein either if the Ricci operator p has only one real eigenvalue 
A or if the Ricci operator p has two complex eigenvalues Ai, A2 with Ai = A2. This does 
not imply that p is diagonalizable in the higher signature setting and hence 971 need not be 
Einstein. 

Theorem 3.3. Let 97t = (V, (•, ■), A) be a 0-model of arbitrary signature. If Wl is Puffini- 
Videv, then we may decompose Wl — Wl\ ® ■ ■ ■ ® Wlk as the direct sum of pseudo- Einstein 
0-models 971;. 
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